Equivalence between XY and dimerized models 
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The spin- 1/2 chain with XY anisotropic coupling in the plane and the XX isotropic dimerized chain are 
shown to be equivalent in the bulk. For finite systems we prove that the equivalence is exact in given parity 
sectors, after taking care of the precise boundary conditions. The proof is given constructively by finding unitary 
transformations that map the models onto each other. Moreover, we considerably generalized our mapping and 
showed that even in case of fully site dependent couplings the XY chain can be mapped onto an XX model. 
This result has potential application in the study of disordered systems. 

PACS numbers: 75.10.Pq, 75.10.Nr, 75.10.Jm 



Exactly solvable models play an important role as limit- 
ing cases of more complex system or for testing numerical 
algorithms. Moreover their physical properties can generally 
be calculated exactly and traced back to simple mechanism 
that can be used in more complicated scenarios. In this article 
we consider two notable solvable models: the anisotropic XY 
model -originally introduced in [ 1] with the aim of gaining in- 
sights on the long range properties of the Heisenberg model- 
and the dimerized XX model -used sometimes as a proto- 
type model to describe spin-Peierls distortion-. We prove 
the equivalence of these two models, despite in the literature 
they are generally considered as separate. The equivalence 
is shown directly by means of a unitary transformation for 
their fermionic counterpart and traced back to the spin mod- 
els, carefully taking care of the boundary conditions. 

Introduction For a chain of length L the dimerized XX 
and anisotropic XY models are given by the following Hamil- 
tonians 

fl3 = 5E(i+7(-i) 4 ) [«+i + ^f+i] a) 

H xy = \ E [(! + T) «+i + C 1 - 7) "Wfi] • ( 2 ) 



The superscript 77 denotes different kind of boundary condi- 
tions (BCs), cr^'i = f]<f\ V , f) = 1,-1,0, corresponding to 
periodic (PBC), antiperiodic (ABC), and open (OBC) respec- 
tively. Since both Hamiltonians Eqs. (Q~|) and (J2J commute 
with the parity operator P = Jli^f' we define the parity 
sectors a = ±1 and the corresponding projection operators 
II ff = (1 + crP) /2. The central result of this article estab- 
lishes that Hd and Hxy are unitarily equivalent (=) up to at 
most a border term as precisely stated by the following 

Theorem 1. For L odd and OBC the models (Q and (0 are 
unitarily equivalent. For L even and PBC or ABC the equiva- 
lence holds in given parity blocks depending on the boundary 
conditions, according to the relation: 
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In other words, for L even, the boundary index in one model 
sets the parity sector in the other (times a modulation factor 

(-l) i/2 ), i.e. a X Y = {-l) L/2 rj d and a d = {-l) L/2 rj X Y- 
An immediate consequence of this result is that the two mod- 
els share the same thermodynamics, since for L — > 00 the 
effect of boundary terms disappear. 

The reason for considering OBC is partly due to the pos- 
sibility of using models Eqs. (Q} and (O to implement quan- 
tum information devices. It has been shown in [2] that in the 
ground state of the dimer model (though with OBC and L 
even), the end spins tend to entangle considerably already for 
small values of the dimerization 7. Moreover the entangle- 
ment survives in the infinite length limit (long distance en- 
tanglement). In a similar fashion, it was already observed in 
III] that the end-spins of the anisotropic model order and such 
order survives in the thermodynamic limit (TDL). However 
this kind of order is of classical nature and no entanglement 
is present between the end-spins of the open anisotropic chain 



Before proceeding to examine the proof of the Theorem, 
let us spend few words on some benefits of such result. First, 
let us note that both models commute with n -rotations around 
axis x and y, H.% = e iTrcr ^ / 2 , a — x, y. However, the 
dimer model Hd manifests a much larger symmetry, the total 
magnetization M z = J2i °f ■ This means that Hd is block 
diagonal in sectors with given magnetization M z , a feature 
which is especially useful in case of non-integrable extensions 
of Hd (which maintain this symmetry) where one has to re- 
sort to numerical diagonalization. Thanks to TheoremQ] such 
a symmetry (or an approximate one) must exist also for the 
anisotropic model Hxy- As we will see the magnetization in 
the dimer model is mapped onto a non-local operator which 
we are able to compute. Clearly this operator has the same 
spectrum of M z and commutes with Hxy- 

The proof of Theorem[TJrelies on a similar theorem holding 
for the fermionic version of the models (denoted here with a 
tilde), 
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Here e = 1, -1, distinguishes among PBC, ABC, and OBC 
for the fermions, i.e. d^ + i = ed% and <Xl+i = eai. As we 
will see later, the spin systems (0 and (0 are connected to 
the quadratic fermionic models and H XY y i a a Jordan- 
Wigner (JW) transformation, after careful reshuffling of the 
boundary conditions. The result for the fermionic models is 

Theorem 2. In the following cases: L even and PBC or ABC, 
L odd and OBC, the models (0 and (0 are unitarily equiva- 
lent, i.e. there exists a unitary operator U (a "mapping") such 
that UH e d W = H XY . 

A simple way to remind the different cases in which the 
theorem applies is given by the following argument. Sending 
aj — > idj in H XY , one realizes that the spectrum of H XY is 
invariant under the transformation 7 — > —7. By relabeling the 
sites of the dimer model, one sees that possess the same 
invariance only when it contains an even number of bonds. 
This occurs for L even in case of PBC or ABC, while for L 
odd only in case of OBC. 

Proof of Theorem 2 Since the fermionic Hamiltonians are 
quadratic, one way of proving the equivalence between them 
is to show that they have the same one-body spectrum. To di- 
agonalize the anisotropic model we rewrite the Hamiltonians 



following the conventions of (Tl]: H 

(V2) feij + h-c.1 and Hf, 

matrices given by 
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The one particle energies of H XY 816 given by the (positive) 
square root of the eigenvalues of (A — B) (A + B). Calling 
Afc such roots, since A is traceless, one arrives at |[Q] 

h X y = ^k-nim - - A fe- 

k k 

The equivalence of the two models now stems from the fact 
that, for L even and PBC or ABC, and for L odd and OBC 
M 2 = (A — B) (A + B). Moreover, under the same hypoth- 
esis, the eigenvalues of M are symmetric around zero (for 
L odd and OBC there is one zero eigenvalue). To write H\ 
in the same form as H XY perform a particle-hole transfor- 
mation on the negative eigenvalues of M. We arrive then at 

H d = Efe - E„og A *' where E„cg is the sum over 

the negative eigenvalues of M. To complete the proof note 

that, in the specified cases, Encg ^ k = (V^) Efc ^-k- D 



The mapping The above proof does not give the explicit 
form of the mapping. We will now provide a physically more 
compelling proof which has the additional advantage of re- 
vealing an exact form of the mapping. For simplicity we will 
stick to L even and PBC/ABC for the fermionic models. The 
first step is to write both models in Fourier space 



H% = [2 cos(fc)<4<ifc + 2ij sin (k)d\. +n dk 

k 

|2 cos{k)a\,a-k 
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Let us consider first PBC. The momenta in the Brillouine zone 
(BZ) are given by k = 2nn/L, n = -L/2 + 1, . . .,L/2. 
Note that only for PBC and ABC if k e BZ then —k e BZ. 
Moreover only for L even k E BZ =>■ k + tt E BZ. In 
particular Eqs. © and © are not correct if L is odd. The 
unitary transformation that maps the dimer model onto the XY 
is 




CL-k-ir < k < 

< k < TT, 



(7) 



Notice that the particle hole transformation does not involve 
neither k = nor k — tt. In fact, for these two momenta, the 
dimer model is given by 2(dJdo — d^-d„) and the anisotropic 
one by 2(a^a — a[a T ). The same mapping Eq. (0 transforms 
H\ into Hx Y a ^ so m ine case °f ABC where the momenta 
satisfy fc = tt/L (2n - 1), n = -L/2 + 1, . . . , L/2. 

The mapping Eq. (jTj can be written in a compact form as 
dt = f + (k) at + /_ (k) a_fe_ 7r with the help of two auxiliary 
functions f± (k) := 6[± sin (k)] ± <S S in(fc),o/2, where 9 is the 
Heaviside function with the convention 9 (0) = 1/2. 

Thanks to Eq. (0 the equivalence between (fermionic) 
dimer and anisotropic models can be generalized. In fact, the 
mapping transforms an r-nearest neighbor hopping term into 
itself, provided r is odd. Instead, an alternating hopping of the 

formE; (-1)* d\d l+r + h.c. becomes Ei( a I a l+r + ai+ra t ), 
again for r odd. When r is even the mapping introduces non- 
analiticities in Fourier space and correspondingly the trans- 
formed model becomes long-ranged in real space. These find- 
ings can also be obtained directly in real space Fourier trans- 
forming back Eq. (0: 



4n = [f+ ( TO ~ x ) a * + (- 1 ) 35 /- ( m _ x ) a i 



, (8) 



with the definition f± (x) = L^ 1 E fc e~ ikx f± (k). Writing 
simply f± in place of the matrix (/±)i,j f (i — j) the fol- 
lowing relations hold: /±/± = /±, and /+/_ =/_/+ = 0. 

Proof of Theorem 1 The first step is to map the spin mod- 
els Eqs. and (0 to fermionic models via the JW trans- 
formation. In terms of ladder operator af = (erf ± ia y ) /2, 

the JW is given by of = cje t7r ^j =1 c j Cj (this in turn implies 
cr,r = Cie" i7r ^=W c J, erf = 2c\ci - 1). The dimer and 
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anisotropic boundary terms become respectively 



introducing Majorana fermions ( a (j), ol = 1, 2, 
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where A^( a ) is the total number operator for the d (a) 
fermions and tjmxy) specifies the spin BC for the dimer 
and XY model. For OBC, rjd = Vxy = 0, we can di- 
rectly apply the result of Theorem [2] and deduce that also 
the spin models are unitarily equivalent for L odd. To study 
the remaining cases we first need to compute exp (inNd) un- 
der the action of the mapping Eq. (0. Writing the number 
operator in Fourier space we get Nd = X)o<fc<7r a \ a k + 

E-7r<fc<o a k a k- Th e sum over negative momenta contains 
a different number of terms depending on the boundary con- 
ditions. For PBC the sum contains L/2 — 1 terms while for 
ABC it contains L/2 terms. Calling N+ = J2o<k<n a k ak 
and N~ = J2_ n<k<0 a l a k we can write compactly N d = 
N+ - N~ + L/2 - (1 + e) /2, where e = ±1 defines the 
boundary conditions of the fermions. Since are integers, 
under the action of the mapping we obtain exp (inNd) — 

L II 

— e(— 1) exp (inN a ). Let us now consider the spin mod- 
els H^ d in the parity sector ad (cxy)- Thanks to 
the JW transformation, for L even, the parity operator is 
P = e l7rAr<1 ( a > and so in each sector e i7rAr<i ( a > = ctmxy)- 
This means, that in the parity sector ad the spin model with 
BC T]d has boundary conditions —r\d<3d in the fermions. The 
same clearly holds for the anisotropic XY model. The equiv- 
alence of the fermionic models (Theorem|2]i holds when they 
have the same BCs, that we denote with e. So, we arrive 
at the relation J]d<?d — Va&a — — e. Now we use the map- 

L II 

ping of exp (inN d ) = cr d , obtaining a d = -e (-1) a X y- 
Solving these last two equations, we finally obtain the par- 
ity sectors and boundary conditions, under which the equiv- 

L 11 

alence of the spin models apply: cfxy = f)d (— 1) and 
Vxy = o d (-1) L/2 . □ 

In the above proof we have partly seen what happens 
to the conserved quantity Nd after the action of the map- 
ping. The precise form also depends on the boundary 
conditions e. In Fourier space we can write Nd — 

E k [/+ ( fc ) - /- M 4 a k + L/2 - (1 + e) /2. Since the 
functions f± (k) are not analytic, and using /+ (fc) = 

1, Vfe, the number operator becomes non local in real space. 
For example its explicit form for PBC (e = 1) is 



2 L v^/^-i/ 2sin[(a;-y)(f + f)l 
N d = j-N a +--1+J2(i) V 1 -'- 

x=jty 



Lsin [(x-y) f] 



~ a l a v 



(9) 

A corollary of our proof is that such operator commutes with 
the anisotropic Hamiltonian Eq. © (e = 1) and its spectrum 
is made of integers from zero to L. 

Majorana fermions The possibility of mapping the dimer 
model into an anisotropic one is not restricted to the mapping 
Eq. (|7]i. Another mapping is obtained directly in real space by 



Ci(i) 
C 2 (i) 



-( 1 1 



that satisfy the commutation relations {Ca(?')> Cpti')} = 
8 a p5jji. This way, it is possible to show that each model gets 
transformed onto two separate Ising chains in transverse field 
each consisting of L/2 sites yfl. Then, assuming L even and 
PBC or ABC, the two pairs of Ising chains are made identi- 
cal by translating by one site one of the two chains obtained 
from Hxy- The composition of all these steps yields to the 
following mapping 



4 = 2 



1 r 



"3+ 



(10) 



The transformation above has the advantage of being local in 
real space and much simpler than Eq. (HJ. By using Eq. ( TTOb 
one can reproduce the results of Theorem[2]for PBC or ABC. 
However Eq. ( fTOb is more powerful in view of its applica- 
tions to more general local Fermi models. Using the mapping 
Eq. ( fTOb the "disordered" tight binding model 



H d = ^ Jjd'jdj+i + h.c. 



(11) 



with arbitrary hopping rate Jj can be mapped onto the gener- 
alized anisotropic model 



H X y = ^ [ J j +S> a j a J+ 1 + J j~M a I+i 



h.c., (12) 



with Jj = (±y(Jj ± Jj_i)/2 and Jj must always be con- 
sidered periodic i.e. Ji+i = J%. This mapping can be further 
generalized by adding a uniform and staggered chemical po- 
tential. After applying the transformation Eq. ( fTOb such terms 
become 



L 

[m + Mst (~1) J d)dj - fi- = 



L 



Yl [^+M St i-iy 



a]a J+ i - (-1) J a]a] +1 



-h.c. 



The equivalence between the generalized models Eqs. (TTTb 
and ( fT2b has potential applications in the study of disordered 
systems. To obtain results on the random version of the 
anisotropic model Eq. (fT2l . it can be favorable to simulate 
Hamiltonian Eq. (TTTb which conserves the number of excita- 
tions. Moreover, through the JW transformation, apart from a 
possible border term depending on the BCs, the equivalence 
between Fermi models can be extended to their spin counter- 
part. In this way a random XY model can be mapped into a 
random XX model. 
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Continuum limit The mappings that we have analyzed so 
far admit a simple interpretation in the continuum limit. To 
this end we expand the fermionic fields into chiral components 



ip(x) = 



a ikp-x 



R(x) 



~—iki 



x L(x). For 7 = the two models 



merge in free massless fermions: Hq = ^ . atoy+i + h.c. 
where the band is half filled, so kp = ir/2. In the continuum 
limit, we get 10, H 

H = i [ dx[:R\x)d x R(x)~ L\x)d x L{x)-] (13) 



while the mass-generating terms in 7, Oxy — a ] a ]+i + h- c - 
and Od = {— lp'ajdj+i + h.c. become 

Oxy = i :L\x)R^{x) - R(x)L(x): 
O d = i :L\x)R(x) - R\x)L{x): . 

From these expressions, we see directly that the terms multi- 
plied by 7 in Hd and Hxy are transformed into each other by 
particle-hole exchange (and a minus sign) on the left movers, 
L — > — L', which is reminiscent of the discrete mapping 
Eq. (0 where the particle hole transformation was also ap- 
plied only for negative momenta. 

Translating into bosonic language, it is known that the 
model Eq. ( fT3l is equivalent to the Gaussian model 

H ° = \J dx{[d x e(x)} 2 + {d x $(x)} 2 } . (14) 

The fields $ and are bosonic and reciprocally dual: d x & = 
d T Q and 9 T $ = d x Q. A nonvanishing value of 7 has the 
effect of transforming Eq. ([Pil l in the sine-Gordon model by 
adding a relevant (in the renormalization group sense) term 
Oxy = :sin(\/47r6(x)): or Od = :sin(v / 47r < I ) (a;)):, respec- 
tively in the XY or in the dimer case. Hence, in the bosonic 
language, the dimer<H>XY mapping simply acts by swapping 
<f> <H- 0. It is interesting to observe that a direct consequence 
of the mapping is the interchange between density and current 
density, as it can be readily inferred by their expressions 

p(x) = :RHx)R(x)+LHx)L(x): = — \=d x $(x) 

j(x) = :R\x)R{x)~ L\x)L{x): = -=d x Q{x). 

Integrating these densities over the space we obtain two quan- 
tum numbers: the total number and the current. In particular, 
the total number and current are directly related to the two 
winding numbers m, n € Z of respectively and $. Such 
integers ( which determine the scaling dimensions of the pri- 
mary operators in the Gaussian model) are both good quantum 
numbers for 7 = 0. For 7 7^ 0, the breaking of translational 
symmetry in the dimer chain invalidates the conservation of 
the current, but maintains the particle number conservation. In 
the XY model the situation is just reversed: the particle num- 
ber is no more conserved, due to the pair creation-destruction 
terms, while the current keeps being a good quantum number. 



Higher dimensions The mapping described in Eq. (0 can 
be easily generalized to /^-dimension. In an hypercubic D- 
dimensional lattice the anisotropic model reads 



Hxy — 



h.c. 



(15) 



where x = (x\, . . . , Xr>) and e, is the unit vector along the 
i-th direction. We do not specify BCs here, to fix ideas we 
can take PBC on a bipartite lattice. After Fourier transform- 
ing one realizes that the BZ is contained in [— it, tt] d . Now, 
let us divide the BZ in two regions according to the sign 
of the first moment fci: A = {k E BZ. : k\ E [0, tt]} and 
B = {k E BZ, : ki E {—tt, 0)}. The canonical transforma- 
tion 4 = 4 for k E A and = rf-k-ff for k E B with 
7? = (71-, 7r, . . . , tt), generalizes the one-dimensional version 
Eq. (0. This mapping transforms the Hamiltonian Eq. (TT~5b 
into the D-dimensional dimer model: 



H 



XY — 



[2cos(k-ei)44c 

i=l k 

+ ^fysin(k • e^d^k+ff + h.c. 

f ^(l+ 7 (-lf l )4dx + e,+h.C., 



where the modulation factor is given by (— = 
exp (in ■ x). 

Conclusions In this paper we have analyzed two common 
spin models (XY and dimerized XX) and showed that they 
are unitary equivalent apart from at most a border term. By 
explicitly providing the unitary transformation we have been 
able to generalize the equivalence in many ways. For exam- 
ple the fully disordered (with site dependent couplings) XY 
chain can be mapped onto a disordered XX chain. Consider- 
ing the fermionic counterpart we have also shown that gener- 
ally a dimerized, r-nearest neighbor, hopping term, is mapped 
onto an r-nearest neighbor pair creation term. In one dimen- 
sion our mappings have a simple interpretation in the contin- 
uum limit in terms of bosonic fields. Similar considerations 
can also be extended to higher dimensions. 

Mapping XY models onto XX ones can be useful in view of 
numerical simulations of disordered models or non-integrable 
extensions. This is due to the explicit particle number con- 
servation of the XX models which makes them easier to treat 
numerically. A by-product of our analysis is that particle num- 
ber symmetry is also present in the XY models although in a 
hidden fashion. 
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